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The commonly made assumption that the distribution of current along a dipole antenna is sinu- 
soidal is examined critically together with the associated distribution of charge and the electromag- 
netic field. Measured distributions of current and charge are compared with those assumed in the 
sinusoidal theory. It is shown that the distributions of charge differ so greatly that it must be con- 
cluded that the electric field near and along the surface of the actual antenna is not well approximated 
by the near field obtained with the sinusoidal assumption. 

A new approximate solution of the integral equation for the current in a center-driven cylindrical 
antenna is obtained in terms of simple trigonometric functions. It is shown that this and the associ- 
ated distribution of charge are in much better agreement with measurements than results obtained from 
the sinusoidal assumption. Approximate expressions are then derived for the electromagnetic field 
near the antenna in a form that satisfies the boundary conditions. 


1. Introduction 


Studies of dipole antennas often postulate the current 
1,(z) =In sin ko(h— |a| (1) 


where h is the half length of a perfectly conducting cylinder of radius a. The axis of the antenna 
coincides with the z-axis of the cylindrical coordinates p, 0, z; ko =27/k=w/c where c= V po€o 
is the velocity of light. The distribution (1) has been used not only for far-field patterns, but to 
determine the near field and the admittance. Recently this procedure was applied to antennas in 
dissipative materials including magneto-ionic media in which the surface of the antenna may ac- 
quire a charge [Ament et al., 1964]. 

It is the purpose of this paper to examine the limitations of the sinusoidally distributed current 
and to derive an alternative simple distribution that is a better approximation. 


2. Sinusoidal Distribution and Its Implications 


The current (1) and the associated charge per unit length, 


; dl, es 
q(z) ae ee cos kyp(h—z), z20 (2) 


are rigorous for a coaxial line that supports only the TEM mode and is terminated at z=A in an 

ideal open circuit, Z =~ (fig. la). An ordinary open end (fig. 1b) is not an ideal open circuit since 

as the end is approached the capacitance per unit length decreases to half the value which obtains 

at distances h—z> b—a where 6 is the inner radius of the outer conductor and a is the radius of 
the inner conductor [King, 1955a]. As a consequence, q(z) increases to double the value given 

by (2) at z=hA and the slope of the current curve departs from the sinusoid in (1). This capaci- 

tive end effect is confined to distances of the order of magnitude of b—a [King, 1955b]. If the 

outer conductor is removed the inner conductor becomes an antenna (fig. 1c) and the “‘end effect” 
extends for at least a quarter wavelength from the end. 
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In a coaxial line in which (1) and (2) are accurate, the electric and magnetic fields are 


I 
E,\p, 2) =5 ae Balp, j= 2 5 


where vp= l/p. The scalar and vector potentials g(p, z) and A(p, z) satisfy the relations 


’ 


Ey(p, 2) =e. Behe =idilon 


It follows that 


¢(p, 2) --|" E,(p, z) \dp = pene Inf 


2 =i 
Alp, =~] Bip, dp =e? In® 


Note that 


- ; 1 0 j 
Ep, 2) = 2002) judd, =z [2121.0] n& 


With (1) and (2) 


Oqz) ,j@,/\_ 
a ee L{z)=0 


so that E,=0 for all values of p. 


(3) 


(4) 


(5) 


6) 


(7) 


(8) 


The distributions (1) and (2) of current and charge in the coaixal line (fig. 1a) cannot be correct 
when the shield is removed to leave the base-driven monopole shown in figure Ic since the boundary 


conditions are not the same. For the coaxial line E; =— V26(z) at p=a, Ez= 
monopole E,=— V6(z) at p=a and a radiation condition must be satisfied at infinity. 
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0 at p=); for the 
(Note 


h 
that the Dirac delta-function 6(z) has the properties 8(z)=0, z ¥ 0; 6(z) =, z=0; &(z)dz= 1.) 
If a sinusoidal current is postulated for the antenna it must be maintained by a mére elaborate 
system of generators than a single delta-function V¢6(z). Actually a continuous distribution of 
sources is required. 

The electromagnetic field maintained by the sinusoidally distributed current (1) is [King, 
1956a] 


Bop, 7 ere (cit + eThrn —2 cos hohe) 


(9a) 
UU mbo (z a5 zth_.. 2z ; 
E,(p, j= dar p ern e thorin wees e—tkoran aren cos khe-ier) (9b) 
—jlm —Fhor in —Shoren . 


where ro= V2? + p?, rin = V(h—z)? + p?, ron= V(A + 2)? + p? and fo = 1207 ohms. These formulas 


have been derived from the vector and scalar potentials, 


A a ie 7c 
Ap, 2) Vas ie 2(z’) 


where r= V(z—2')? +p? and where /,(z) and q(z) are in (1) and (2), by means of the relations 


dz' (10) 


en-ikr |, nh 1) € hor 
- dz', Ap, z)= q(z eae 


aE 


— 0A: > 
Bip, = aP Biss 


=e: 2) Ep, =e. 2 —jwAAp, z). (11) 


At a small radial distance p given by axp <h, koa <kop <1, the approximations ri, =h—z, 
ren = h+z, ro =z may be made in the exponents. It follows that 


T(z) 


Bo (p,z) = 5 lar Top 


In ap sin ko(h—|z|) = (12) 


Thus, near the: antenna, Bo(p, z) is proportional to the current as in (3) for the coaxial line. 
If short distances near z=0,+h are excluded, and the same conditions are imposed on p, 


q(z) 
27re9p 


=! Tinto 
amp 


E,(p, z)= cos kh —z)=s—— (13) 


where the upper sign is for z <0, the lower for z<0. It is seen that near the antenna E,(p, z) 
is proportional to g(z) as in (3) for the coaxial line. 

Since at p=a the surface densities of current and charge are K,(z)=I,(z)/27a and n(z)= 
q(z)/2zra, it follows that at p=a (12) and (13) reduce to the boundary condition for the tangential 
and the normal electric fields. 

Although Bo(p, z) and E,(p, z) are the same near and on the antenna as on the coaxial line when 
the currents in both are sinusoidal, this is not true of Ep, z). On the coaxial line Ea, z)—V<8(2); 
on the antenna with an assumed sinusoidal current (9c) applies. However, it does not reduce to a 
simple form at p=a. Its magnitude at p=a for antennas with koh == and a and ==7.022x10-3 
is shown in figures 2a and 2b. Evidently, the distribution of E(a, z) obtained from (9c) must be 
maintained along the antenna by a distribution of sources. Near the antenna where p << h and 
kop << 1, the general relation V X E=—jwB gives 


dEAp, 2) _ 6E,(p, z) 2) 1 [= 1 dq) jo 


8z > ig EMP ona es Be +210] =o. (14) 
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Ficure 2. (a) Magnitude of E,(p, z) at p=a for an an- 
tenna with koh = 77/2. 
(b) Magnitude of E,(p,z) at p=a for an an- 
tenna with khh=7. 


This follows when (12) and (13) are substituted for J,(z) and q(z). Near the antenna, where (12) and 
(13) are acceptable approximations, Ep, z) is practically independent of p and equal to the value 
atp=a. Thatis, Ep, z) = EXa, z). 

A continuous distribution of sources along the antenna and the associated nonvanishing 
Ea, z) require an outward radial component of the complex Poynting vector 


1 
Sola, 2)=— 5 vol (a, 2)Bg(a, 2) (15a) 
instead of the value 
1 1 
S,(a, z)= a) vo Aa, z)Bx(a, 2) =5 voV SBR(a, z)5(z) (15b) 


which obtains on a perfectly conducting center-driven antenna. Note that with (3) 


h 
P=Re 2ma | Spa, 2dz=Re 5 Vsf"(0) (16a) 
—h 


which is the total power supplied to and radiated from the antenna with a delta-function generator 
at its center. The total power radiated from the antenna with an assumed sinusoidally distributed 
current is 


A h h 
P=Re 2na { S,(a, z)dz = — Re 27a [ 5 YE a, 2)BF(a, z)dz =— Re +f Ea, z)I*(z)dz — (16b) 
-h -h -h 


where E,(a, z) and Bo(a, z) are given by (9c) and (9a). 
432 


~—— 


S 
AY 
\ 
ry 
‘ 
; 
| 


—— SINUSOIDAL THEORY 

20} © MEASURED (MORITA) a { : 

” ‘ 

mae + 

= 10.12 yf { 
5b kh = % 7 ' 4 

t 

e 

H 
10 { 4 

SINUSOIDAL THEORY HF 

4 
e-e-6 MEASURED (MORITA) 06 on | 
A 
0 Hl eabat oars 1 rt Het 
0 5 1.0 15 50° -100° ~-50°) = 0 
co) 10 20 
me £002 ReLaTwE scaLe (2) 
—— MILLIAMPS / VOLT sf 
Vo 
Ficure 3. Current along a cylindrical half-wave mono- Figure 4. Distributions of charges along a half-wave 
pole (koh=7r/2, h/A=1/4) over an image plane; O=10.12. dipole. 


3. Comparison With Experiment 


In figure 3 the current along a half-wave monopole is shown when 0 =2 In oh 10.12 are = 78. 
a 


The corresponding curves for the cosinusoidal distribution (solid line) are in the form 


TA2)__ cos kez__ cos kuz 
Ve" Der LE a0 


where Zo= 36.56 + 721.25 = 43.3 exp(j30.2°) 2. The general agreement between the measured 
values and the sinusoidal approximation is fair, with more current near the top of the actual antenna 
than is indicated by the cosine curve. Also the phase differs somewhat from the constant required 
by the sinusoidal theory. The large difference between the sinusoidal approximation and measure- 
ment in the distribution of charge along a half-wave dipole is shown in figure 4 where the measured 
amplitude and phase are in broken lines. The equation of continuity, dl, (z)/dz+jwq(z) =0, when 
applied to (17) gives 


cg(z) _ sin koz 
Ve (Zo 


e —i(00+90°). 


(18) 


The solid curves in figure 4 are for |cq(z)/V%| and 0,=— (0+ 90°) =— 120.2°. The agreement be- 
tween them and the measured amplitude and relative phase angle is quite poor. A large part 
of the charges actually present near the driving-point and the end of the antenna is not included 
in the cosine curve. The fact that the measured phase changes rapidly and greatly at some dis- 
tance from the driving point is ignored by the sinusoidal theory. 

Corresponding graphs for driven full-wave antennas are in figures 5 and 6. Note that the 
sinusoidal theory fails completely near the driving point where instead of I,(0)=0, J{0) is actually 
nearly three-quarters of the maximum value. There is also more current in the upper half of the 
antenna than in the sinusoidal distribution and the phase angle, instead of independent of z changes 
significantly along the entire antenna. The maximum current Jm in figure 5 has been chosen so 
that the total powers radiated by the sinusoidal and the actual currents are equal. Specifically, 
P/V? = Go = |Lo/Vo|2Ro = 1.023 X 10-3 mhos is the power per unit voltage squared supplied to the 
actual antenna. The power radiated by the sinusoidal current is 2, R¢ = 199.12. This is equal to 
1.023 x 10-3 when In =2.26 X 10-3V¢. 
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FicurE 5. Distributions ‘of current and charge along Iq(z)} RELATIVE SCALE 
center-driven full-wave antenna; koh=77, O= 9.92. 


FicureE 6. Distributions of current and charge along 
center-driven full-wave antenna; khh = 7, Q=11.5. 


The distribution of charge along the full-wave dipole (fig. 6) varies almost linearly with distance 
so that large concentrations of charge near the driving point and at the ends are not included in a 
cosine curve. The abrupt reversal in phase at z=0.25A, which is predicted by the sinusoidal 
theory, is replaced by a rapid reversal near z = 0.3A. 

In view of the reasonable agreement between the sinusoidal approximation and actual dis- 
tributions in some respects and the large differences in others, the question arises as to how ac- 
curately such quantities as the far and near field, the admittance and the radiated power can be 
determined from the sinusoidal assumption. 


In the radiation field (kor > > 1) the electric and magnetic vectors differ only by a constant co- 
efficient since E4,=cB%. Moreover, both depend directly on the integral 


h 
i T(z’ ete’ £8 © gin © dz’ (19) 
which is insensitive to the details of the distribution of current and is independent of the charges. 
Hence the far field patterns of the actual currents (figs. 3 and 5) should be in reasonable agreement 


F : : nN 
with the patterns of the sinusoidal currents. (In the case of longer antennas (h > 3) the patterns 


involve minor maxima and minima that depend on rather precise relations between the phases and 
amplitudes of the currents along the antenna, so that greater divergences between the measure- 
ments and the sinusoidal theory may be anticipated.) 

When the currents are normalized and the antennas and the surrounding medium are lossless, 
the total power radiated is the same for both distributions. The real part of the integral of the 
normal component of S = }19E X B* over any surface which completely encloses the generator is, 
therefore, also the same for the two distributions. 

Whereas the magnitude of the magnetic field of the sinusoidally distributed current is a fair 
approximation of the magnetic field of an actual cylindrical antenna at all points (since the sinusoidal 
distribution is a reasonable approximation of the current), this is true of the electric field only in 
the far zone. The electric field near a sinusoidally distributed current differs greatly from that of 
a center-driven cylindrical antenna since FE, near the antenna depends on gq(z) and it has been shown 
that the cosinusoidal approximation for q(z) is in poor agreement with the actual distribution. 
Moreover, since a sinusoidally distributed current cannot be generated by a single source along a 


perfectly conducting cylinder, a continuous distribution of generators is required to maintain an 
impressed tangential electric field along the entire antenna. A resultant nonvanishing tangential 


electric field then obtains and this is not a good approximation of F, = 0 along a perfectly conducting 
cylinder. 
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It is not actually possible to define a driving-point impedance for an antenna with a sinusoidally 
distributed current since there is no single localized generator. A quantity that is dimensionally 
an impedance may be introduced in terms of the total power defined in (16b) and the current at 
the center of the antenna. Thus, 


Z- 2P 

° T{0)F*(0) (20) 
Since the generator is distributed along the entire antenna (20) is not a driving-point impedance. 
The driving-point impedance of the cylindrical antenna is 


Zo =VE/TA0). (21) 


It is not surprising that (20) for the sinusoidal theory and (21) for the center-driven cylindrical 
antenna differ greatly except in special cases. Zo in (20) is a function of koh only, whereas in (21) 
Zo depends on both koh and”. When koh =9> (20) gives Zy)= 73.14 42.5 O or Yo=(10.22 — 75.94) 


x 10-3 mhos, whereas (21) ranges from Zo= 78.5 + j43.6 when ne 11,000 to Z)>= 108 +36 when ‘ 


~ 16.5. When kyoh=7, (20) gives Z>=% and Yo=0, (21) gives Zp = 2674 — j2630 when ae 11,000 
and Zo = 264—j436 when g ~ 75. Actually, (20) with a sinusoidal current is a reasonable approxi- 
mation of (21) for a highly conducting cylinder only when koh is near 7 and q = 10,000. Since the 
driving-point impedance depends critically on both the electric and the magnetic fields at the sur- 
face of the antenna and since the sinusoidal theory provides a fair approximation of the magnetic 
field but not of the electric field, general agreement cannot be expected. 


4. Integral Equation and Its Solution 


The measured distribution of current, the impedance, and especially the distribution of charge 
and the near field of a center-driven cylindrical antenna differ significantly from the quantities 
obtained when a sinusoidally distributed current is assumed. More accurate solutions of the in- 
tegral equation for the current have been derived by iteration [Hallén, 1938; King and Harrison, 
1943; King 1956b], by variational methods [King, 1956b; Storer, 1951; Tai, 1950] and by judicious 
separation and fitting [King,-1959]. It is the purpose of this study to derive a new solution which 
combines considerable accuracy with reasonable simplicity for cylindrical antennas with kgh<577/4. 

The differential equation for the vector potential A,{z) at p=a on a cylindrical antenna center 
driven by a delta-function generator is [King, 1956b] 


i 2 
(E+ ka)dto=—-78 Hoe) (22) 


where the time dependence is e, It is assumed that kpa<1, a<h. A-A(z)=Ap(z)+JjAdz) is 
given approximately by 


ree oe us , ’ ’ 
Az) = Ga: [ ; LAz)K(z, 2')dz (23) 


where 


K(z, z')=e er |r, r=V(z—2'? +a?. (24) 
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The real and imaginary parts of the vector potential are 


kyo fh cos kor 
= , t 25 
Anta) = Gat | Bet) EE (25a) 
ies ko is Fi sin kor 7 
A(z) = ives T{z') ae dz’. (25b) 


The solution of (22) is readily combined with (24) to obtain the integral equation for the current 
[King, 1956b]: 


[’ L{2z')K(z, 2')dz' 2 l¢ cos kez + 5 Ve sin hell (26) 
-h 0 


where fo = V po/€o= 1/V vo€o = 1207 ohms. 

Approximate formal solutions of (26) have been described [King, 1956b] and accurate results 
for particular cases can be obtained by numerical methods. In order to derive an analytically 
simple and reasonably accurate solution for [,(z) use may be made of the quite different properties 
of the kernels in (25a, b). In particular, the kernel in (25a) is very large at z’=z where rma, the 
kernel in (25b) is not. Specifically, at z'=z 


cos kor _cos koa . 1 
kor eat (27a) 


sin kor sin koa . 


kor koa ! : (27b) 


Graphs of the two kernels in figure 7 show that Ar(z) must vary with z like Iz). However, since 


sin kor cos kor 


FicurE 7. The functions ie be 


and cos (kor/2). 
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I{h)=0 whereas Aj(h) is small but finite, the difference Ap(z)—Ar(h) should be approximately 
proportional to [{z). That is 


4arvo[Ar(z)-Aath) ]=V(2)1(z)=V ar) (28) 


where Wap is the approximately constant value of V(z) defined at a suitable reference value of z. 

The graph of (sin kor)/kor in figure 7 shows that contributions to A;(z) come almost equally from 
elements of current over a large range of |z’—z|._ Hence, an approximation of A,(z) can be obtained 
from 


sin kor _2 sin kor cos $kor 


a. 1 
kor kor TEED kor (29) 


over the range kor<z as shown in figure 7. With (29), (25b) may be expressed as follows: 


ko fr 1 1 
A(z) = =| T{2') [cos 3 kor1 + cos 9 ker] dz’ (30) 


where m=V (z-2’)*4a?=|z—z'|, m=V (ztz')*+a2=|z+2'|._ The trigonometric functions can be com- 
bined to give 


< —2k I h t 1 , , 
A(z) = Tavs 3 ke | IAz') cos 3 koz' dz’. (31a) 
It follows that since 
A (z)=A(0) cos 3 koz (31b) 


where A,(0) is a constant, A,(z) is independent of the distribution of current. Also 
A(z)—A(h)=A(0)F (32a) 
where 


Fi= cos $koz— cos $koh. (32b) 


An integral equation that involves the vector potential difference [King, 1959] is obtained from 
a rearrangement of (26). It is 


er Not — — Ja 5 e | 

ic Ida! Ke, 2'\de! =? | 5 Vee + UF (33) 

where 
Moz=sin ko(h— |z|) (34) 
Foz=cos koz—cos koh (35) 
Kaz, z') = K(z, 2')— K(A, z') (36) 

and 
i. Agee I " Le K(h, 2!d2’ (37) 
k ko 2 Aer a z ? . 
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The properties, (28) and (32), of the integral on the left in (33) suggest the distribution of 
current 


T(z) = Iy(z) + Iu{z) + Inf) (38) 
where 
Tz) = IvMoz, Ip{z) = IuF oz, Ip(z) = IF, (39) 
and Iy, Iy, and Ip are undetermined constant coefficients. 


If (38) is substituted in the integral in (33) and use is made of (28) and (32) the following relations 
are obtained: 


Iv(e")Kan(z, 2’)d2" = Iy(2)Var (40a) 
[- Tole" )Kar(z, z')dz! + Iz) aun (40b) 
[- Ip(z')Kalz, z')dz' = Ip(z)[Vaor +j¥an1] = In(2)V av (40c) 
il Iv(e')Kal2, z')dz" = (Iv[[)Io2) Var (40d) 
[Flee Kare, 24a" = UllooleV ae (406) 


Similarly, if (38) is substituted in (37) 


[5 tee Kh, 2242! = LV) + Le olh) + 10 Vl. (a1) 


Explicit formulas for the several V functions are in appendix A. 

With (38)—(41) equation (33) is reduced to an algebraic equation in Moz, Foz and F,,. However, 
since approximations have been made, it is preferable to make the substitution in the differential 
equation (19) in the following equivalent form: 


Arie (f + is) [Aclz) — A-(h)] = —titks V¢8(2) — 4arvok2A.(h). (42) 


2 
With (24a, b) and (38)—(40e) it follows that: 

4arvo[A2(z) — AAh)] = IvVarMoz + IuVaurF oz + GIVVart jluVavi + InVav)F4,. (43) 
Differentiation and formation of the left side in (42) results in 


[Left side of (42)] =— 2IyVarko cos koh&(z)—IuVaurk2 cos koh 
: . ] 3 ] 
(itVart HeVat L >Ven) Re (cos 2 koh — 40°85 kez) (44a) 


With (41), 


[Right side of (42)] i Ve8(z) — kU vy v(h) + TuVol(h) + [pV (h)). (44b) 
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When (44a) and (44b) are equated, the coefficients of 5(z) must be equal. The resulting equation 


defines Iy as follows: 


l= janVé 
* boar cos koh’ 


(45) 


The remaining equation reveals the approximate nature of the solution since the left side is a 
function of z, the right side a constant. However, within the restriction koh < 57/4, the slowly 
varying function cos $koz can contribute only a small error in the distribution of current if the 
equation is satisfied at z=+h andz=0. In this way two equations are obtained that determine 


Ty and Ip: 
TyCy t+ In€Cp =1vCy 


IvEu + IpEp = IvEy 


where 
3 1 
Cu=— Vaur cos kyhh +jWVaur G- cos 5 keh) + Vy(h) 
3 1 
Cpo= Van (5- cos 2 koh) + Wp(h) 
3 1 
Cy=— [a G — cos 5 koh) + Yolh)| 
wl 1 
Ey=—Vaur cos koh —j 4 Wav cos 3 koh + Vuh) 
1 1 
Ep = “4 Wan cos > koh + WVp(h) 


ar ee 5 koh —Wo(h). 
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The equations (46a, b) are readily solved in the form Jy =Tyly, In>=Toly, where 


Se eee a 
Ty= CyEp—CpEu QO [VaVv(h)—jVarV(h)] 
Cuky — CvE. : 
Tp= G, E, =z C, E, =—jQ7'{Val[ Var cos koh — Vy(h)] + Vaur¥v(A)} 


with 


Q = Wap [Vaur cos koh = WVr(A)] +jVavV (A). 


5. Distribution of Current and Admittance 


(46a) 


(46b) 


(47a) 


(47b) 


(47c) 


(47d) 


(47e) 


(47) 


(48a) 


(48b) 


(48c) 


With the amplitude coefficients Jy, Iy, and Ip determined the approximate distribution of 


current is 
___JénVs if. = z Ly eer 
[{2)= 5oVan cos koh {sin ko(h —z)+ Tu(cos kpz— cos koh) + Tp(cos 5 koz— cos 5 koh) (49) 
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forz 20. Note that /,—z)=I,z). This formula differs by the term in Tp from an earlier approxi- 
mation [King, 1959]. It more accurately represents the distribution of current especially at some 
distance from the driving point. (Note, however, that the added term is strictly a self-term that 
in no way alters the contributions to the shifted cosine term by coupled antennas in arrays.) The 
admittance of the antenna is 
jar 


Yo ~ SyV ar cos koh 


{sin kyh + Ty(1—cos koh) + Tp(1 — cos 5 hah) (50) 
The distribution of charge is 


27regh c is 1 > 1] 
2 = eT S0l'p as cas 4+ =~ 
q(z) Sancas hah {cos koh—z)+Ty sin koz+ 5 Tp sin 5 her} > 0. (51) 


Note that g(— 2)=— (2). 
When koh is near oe (49) and (50) become indeterminate and alternative forms are more con- 


venient. By a simple rearrangement (49) becomes 


Lz y= ra Ta sin koz— sin koh) -Ti(cos koz— cos koh)— Ty(cos 5 koa — cos poh =0 (52a) 
where 
» — Iu+sin koh fo, Tp r 
tg cos kyh ” > cos koh (52b) 


Explicit formulas are in appendix B. These show that 7;, and Tj, are finite when koh =F The 


admittance corresponding to (52a) is 


Yo= 5 He =| sin koh + T,,(1 —cos koh) —T,(1 —cos ish)| (53) 
The charge distribution is 
2 e : 1 Pe | 
qe)= es {eos koz— Ty, sin kets T;, sin 3 kes} ,z20. (54) 


When the antenna is electrically short so that koh < 1, 


142) = 208 baa(1-F) +5 geT(1-F) |; z=0 (55) 


gy ees aati ri | 


1+ Thos; z=0 (56) 


where T=Ty4+ Tp/4. 
Numerical computations have been made for an antenna with ~=7.022 x 10-* for which 


extensive measurements are available [Mack, 1964]. When kyh=7/2 with 0=8.58 and khh=a7 
with 0=9.92 the following quantities apply: 


koh=F 5 thar = 6.218, Ty = 3.085 + 73.581, T, = 1.061+ 0.025 (57a) 


koh=7r, thar = 5.737, Ty =— 0.117 + 0.114, Tp =— 0.106 + j0. 108 (58a) 
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’ : : 7 
The corresponding currents, admittances, and charges per unit length are, for koh = ‘ 


T(z) 


Ve =|9. 597 cos koz — 0.067 cos = 5 hue +-0.047 — j 2.680 (sin koz— 1) + 8.269 cos koz 
0 


— 2.843 cos shoe + 2.010|| x 10-3 (58a) 


Yo= (9.577 —j4.756) X 10-3, Zy = 83.76 + 41.60 (58b) 
m2) {2. 680 cos koz — 8.269 sin koz + 1.422 sin = 5 hoz — ilo. 597 sin koz 
0 
— 0.033 sin we || x 107%. (58c) 


For koh = 77 they are 


a {o. 331 (cos koz+ 1) + 0.314 cos © hye — ile. 905 sin koz — 0.340 (cos koz + 1) 


— 0.308 cos she} x 10-3 (59a) 


Yo= (0.976 + j0.988) x 10-3, Zo = 506.0 —j512.2 (59b) 


cate) =—{[0.331 sin kz + 0.157 sin td 43] 2.905 cos koz +0.340 sin koz 
0 


+ 0.154 sin al x 1073. (59c) 


Graphs of [2 (z)/V§= 19 (z) +jL,(z)]/V6§ are in figures 8 and 10 together with the measured results 
of Mack [1959]. Corresponding curves of cq(z)/Vé= Ge/V5) [q"(z) +jq'(z)] are in figures 9 and 
11. These were obtained from the slopes of the measured curves of the currents. The agreement 
between the approximate theory and experiment both for currents and charges indicates significant 
improvement over the sinusoidal theory. 


TT 
E(z)_ Iz(z)_1,(z) 
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[ i= 2fn "= 854 ‘\ === FROM sore oF measuneD cunves | Vor ale V4 
= r N FOR Att 720 kyh = wh 4 7 
5 = 7022x107? te Tlaz)=T"(2}+jT'tz) 2 4 
h=% ‘ a=2tn 2h 954 a, 
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FIGURE 9. cq(z)/Vé for an antenna with kgh=a/2, O=8.54. 


Ficure 8. 1,(z)/Vé¢ for an antenna with koh = 7/2, O=8.54. 
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FicurE 10. 1,(z)/V§ for an antenna with koh=7, O=9.92. 


6. Electromagnetic Field 


If the three-term approximations of the currents and charges are substituted in (10a, b) and 
these are then used in (11) expressions are obtained for the electromagnetic field that are better 
approximations than (Qa—c). Unfortunately they cannot be expressed in as simple a form. How- 
ever, an approximate expression for F, near the antenna may be obtained from 


I, 
Bap, z) = ane E,(p, z) = we (60) 


which are valid for all distributions of current and charge. It follows from (14a) that 


dEAp,z). 1 [1 qe) “1.0)} 
ap np le az Mm T2) (61) 
With (49) and (51) and for z = 0, this leads to 
= tag |? koh —T Di sas koh)| In®- 
Exp, =P COmtdh uv Cos ko D = cos 5 0z cos 5 0 (62) 
Similarly, when koh is near 77/2, 
— he 
EAp, 2)= Van | sin koh +T;,, cos koh +T5(3 cos = 5 haz cos 5 koh) In . (63) 


Note that these expressions satisfy the boundary condition E,=0 at p=a. They are useful where 
the approximations, r1= V(h—z)?+p?= |h—z|,r2x=VAtz)?+p?=h+z, r= V22+p?= |z| 


are acceptable. 
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For the two antennas to which (57a, b) apply 


koh= z Mp2 A a — | 0.253—j0.018 + (0.809+0.019) cos hoz [in (64a) 
koh =r EsP2 2 fp.) [o. 145 — j0.138 + (0.102 — jo. 095)c0s 5 kee | in®. (64b) 
0 
For comparison with the sinusoidal theory it i x onvenient to use the quantity |AE.(p, 0)/Zolm|. 
With koh=5, 1 = |1,(0)|; with koh = 77, Im = [Iz I The relevant quantities are 
AE Ap, 0) _ AE Ap, 0) Zo _ : p 
kh =", 0.234 +j0.117) In& 
ea UO. gc ieee eng G2) 
nto.3) r6(0.2) v 
koh =a, <> =— =(—0.154+ 0.269) In 2. 
m= TT) OV Cpe a ae (66) 


These formulas are useful when (9/h)? << 1 or (p/A)? << 0.25, which is precisely the region where 
E, as evaluated from the sinusoidal theory is incorrect. 

A graphical comparison of |AE2(p, 0)/Col(0)| for the sinusoidal theory as determined from (9c) 
with the three-term formula (65) is shown in figure 12. The broken line, representing (65), is a fair 
approximation when a is in the range from 1 to about 10 (which corresponds to > in the range from 


0.007 to 0.07 to satisfy the condition p?/A?=0.005 <0.25). On the other hand, at larger distances 
(9 ~ *) the fields given by the three-term theory and the sinusoidal theory should be quite compa- 


rable if I,(0) =Im is the same. Hence, the solid-line curve in figure 12 should be a good approxi- 
mation when an 50 (corresponding to ta 0.35). Between these ranges neither formula is cor- 
rect, but it may be assumed that they join continuously and smoothly as suggested by the dotted 


curve. 
Until an accurate evaluation of the near field of a center-driven antenna is available, a reason- 


able approximation of E,(p, z) is given by (62) or (63) when © is small and by (9c) when r is large. 


In an intermediate range an estimated smooth connection between the two theories may be made in 
the manner suggested in figure 12. 


wane [XEL| oo 76 tn (A/a) (THREE-TERW THEORY) 
Jotz(o) 
AE, AE, (p,0) 0) 


%13(0) Fors keep? 


kyh = 
4 = 7.082 x03 
zrzQ 


(SINUSOIDAL THEORY) 


FIGURE 12. E,(p, 0) near a half-wave dipole. 
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7. Conclusion 


The assumption that the distribution of current along a center-driven highly conducting 
antenna is sinusoidal is inadequate to determine the distribution of charge on and the electric 
field near the antenna. A new relatively simple representation of the current and associated 
charge has been determined which is a much better approximation and which provides a means 
for determining the approximate near field. The: results suggest that approximate solutions of 
dipole antenna problems that depend on an assumed sinusoidal distribution of current and involve 
the environment near and on the surface of a highly conducting antenna may be of questionable 


validity. 


Note that (62) and (63) do not take into account of contributions to E,(p, z) by Coulomb-type fields 
near the generator at z=0 and near the ends atz=+h. These latter are responsible for the large 
peaks in figures 2a and 2b. 


8. Appendix A 
The several VY functions introduced in (40a—e) are defined as follows: 
War(0), koh < w/2 
War = » (A-1 a) 
Warn), 2 heh 
where 
We 11 [cos kor cos korn 1 
Wanlz)=csc kelh—|z|) | sin ko(h—|2"l) [cos-Fer— c08 For dz (A-Ib) 
-h r Th 
Waur = Vavr(0) (A-2a) 
where 
7 cos kor cos kor; ; 
Wavr(z)=([cos koz—cos koh]! | [cos koz’ — cos koh] [sos Aor 208 fora dz (A-2b) 
-h 
Pap = Van(0) (A-3a) 
where 
-1 fh —jkor a Le 
Wan(z) = | cos 3 koz — cos Z lah] i [cos 1 koz’ — cos a koh E ages | dz' (A-3b) 
2 2 —h 2 2 Th 
War = Va(0) (A-4a) 
where 
-1fh ‘ : 
WValz)=— [eos ; koz — cos 3 koh | sin ko(h— |z’|) [ees dz’ (A-4b) 
—h r Th 
Wau = Vaui(0) (A—5a) 
where 


-1 fh P ; 
Wavilz)=— cos s koz — cos 5 koh | , [cos koz’ —cos koh] | seen oe dz’. (A-5b) 


Th 
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A graphical representation of the several functions V(z) for koh = Sn and 7 indicates that these 


functions are in fact reasonably constant at the indicated values z. 
The functions V(h) introduced in (41) are defined as follows: 


e thorn 


dz' 


h 

winy= |" sin ko(h — |z'|) —— (A-6) 

Wh) = [ (cos koz' — cos koh) ii dz' (A-7) 
h ] e—thorn ; 

Wp(h) = ie (cos 5 > hoz’ —cos 5 La dz'. (A-8) 


In the above formulas, 
r=V(z—2' +a’, rh= V(h—2z'? + a2. (A-9) 


All of the above integrals can be expressed in terms of the tabulated generalized sine and cosine 
integrals [Staff of Computation Laboratory, 1949] without difficulty except the following one which 
involves some manipulation: 


h r =e 

i wenn ea — dz’ =| sinh 2 aapeee ‘| ae 

-h 2 a a 2 

1. 1 1 a8 
a5 ella Riu aene 5 halh-+2)) +8 (4' 5 lagh—2)) +35 (4', 3 halh-+2)) | 

ERE ie ores 3 teal £3 
— Fete 1C (4 , 5 hoth-+2) +0 (4 , 5 kah—2)) +38 (4 , 5 hoth+2)) +58 (4 Si) 
(A-10) 


where A’ = (koa V3)/2, C and S are tabulated generalized cosine and sine integral functions. They 
are defined as follows: 


Sla, x)= i van du (A-11) 
C(a, x)= i deco (A-12) 
where W = Vu? +a’. 
9. Appendix B 


The coefficient Tj, is defined in (52b) with (48a). It may be expressed as follows: 


WVaVy(h) —jVarVo(h) + Q sin koh 


T= _ 
u Q cos koh (B-1) 


where 
O=WVan[Vaur cos koh —Vo(h)] +7Vaur P(A). (B-2) 
It is readily seen from (A—1), (A—2), (A-4), (A-5), (A-6), and (A-7) that when koh =» Waur= War, 
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Wau = Var, Vo(h)= Vuh). It follows that 


oT Vato m4 (B-3) 


Similarly, 7; as defined in (63) with (47b) is given by 
; _t{Var [Waur cos koh — Vuth)] + VaurVv(A)} 


a Q cos koh (B-4) 
When kyh= > this reduces to 
_ IVarVaur , LX 
ip= Wan¥ y(h) —jVar¥ a(h) h 4 (B-5) 
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